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Abstract 

We reconstruct B° — > D~D + decays using a sample of 275 x 10 6 BB pairs recorded by the Belle 
experiment, and we measure the branching fraction to be [7.42 ± 0.23(stat.) ± 1.36(syst.)] x 10~ 3 . 
We also search for the related decay B° — > DfD~. We observe no statistically significant signal 
and set an upper limit on the branching fraction of 2.0 x 10~ 4 at 90% C.L. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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Introduction 

Several decay modes of B mesons with in the final state have been measured at 
the B factories. The amplitudes governing these decays are interesting because neither 
constituent flavor of the is present in the initial state. For example, the recently-observed 
decays B° — > DfK~ [1] and B° — > Dfj(2317)K~ [2], with branching fractions in the range 
1CT 5 — 10 -4 , can proceed via a bd — > cu H^-exchange diagram. Here we study the related 
decays -B — > DfD~ and B° — > D~D + . The former proceeds via Cabibbo-suppressed W- 
exchange and has not yet been observed; a recent theoretical calculation predicts a branching 
fraction of ~2.5 x 10~ 4 [3]. The latter proceeds via a Cabibbo-favored tree diagram; the 
ratio of its branching fraction to that for B° —> DfiT~ can be used to determine the ratio of 
CKM matrix elements \V ub /V cb \ [4] assuming factorization [5]. However, the PDG value [6] 
for B(B° — > D+D~) has a large uncertainty (38%), which limits the usefulness of this method 
at present. 

In this paper we report an improved measurement of B° — > D~D + decays and a search 
for B° — > D+D~ decays with the Belle detector [7] at the KEKB asymmetric energy e + e~ 
collider [8]. The results are based on a 253 fb _1 data sample collected at the center-of-mass 
(CM) energy of the T(AS) resonance. To study backgrounds, we used a Monte Carlo (MC) 
simulated sample of T(4S) — > BB events and e + e~ — > qq (q = u, d, s and c quarks) contin- 
uum events; the sample size was about twice that of the data. To evaluate the reconstruction 
efficiency, we used a large MC sample of several thousand signal events. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a multi- 
layer silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scin- 
tillation counters (TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) 
crystals located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. 
An iron flux-return located outside of the coil is instrumented to detect K\ mesons and 
to identify muons (KLM). The detector is described in detail elsewhere [7]. Two different 
inner detector configurations were used. For the first 152 million BB pairs, a 2.0 cm radius 
beampipe and a 3-layer silicon vertex detector were used; for the latter 123 million BB pairs, 
a 1.5 cm radius beampipe, a 4-layer silicon detector and a small-cell inner drift chamber were 
used [10]. 

Charged tracks were selected with loose requirements on the impact parameter relative 
to the interaction point (IP). We also required the transverse momentum of the tracks to 
be greater than 0.2 GeV/c in order to reduce low momentum combinatorial background. 
Based on the optimization of the product of the efficiency and sample purity as evaluated 
by MC, we impose the following requirements. For charged particle identification (PID) 
we combined information from the CDC, TOF and ACC counters into a likelihood ratio 
£(i^ ± )/(£(K ± ) + £(71^)) [11], which we required to be larger(smaller) than 0.5 for charged 
kaon(pion) candidates. This requirement selected kaons and pions with average efficien- 
cies of 92% and 96%, respectively. Neutral kaons (Kg) were reconstructed using pairs of 
oppositely-charged tracks that have an invariant mass within 30 MeV/c 2 of the nominal K° 
mass, and a displaced vertex from the IP. We used the — > <pn~ , K*°K~, and K S K~ 
modes to reconstruct Dj mesons and D + — > K + K~n + , K~it + it + , and K s n + for the D + 
mesons [12]. Combinations of oppositely-charged kaons with — M K + K - \ < 20 MeV/c 2 
and of oppositely-charged kaons and pions with \Mx*o — Mx+ n - 1 < 85 MeV/c 2 were retained 
as and K*° candidates, respectively, where and M K *o are the nominal masses of the 
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FIG. 1: AE (left) and M bc (right) distributions for reconstructed B° — > D~D + events in the 
M bc or AE signal region, respectively. The curve shows the results of the fit. The normalized 
distributions for the events in the sideband of D s and D meson masses are shown as the hatched 
histograms. 

two mesons [6]. These combinations, as well as K^s, were then combined with a or ?t ± 
to form Df mesons. All combinations with invariant masses within a 4 <r(4.5 a) interval 
around the nominal Df(D ± ) mass were considered for further analysis. The values of a 
were determined from invariant mass distributions of MC Df{D ± ) decays and were in the 
range 3.6 - 4.2 MeV/c 2 . 

Pairs of Dj and meson candidates were combined to form B° meson candidates. 
These were identified by their CM energy difference, AE = E^ M — E™ m , and the beam 
constrained mass, M bc = \J {E™ m ) 2 — (Pb M ) 2 , where E™ m = y/s/2 is the CM beam energy 
and E% M and p§ M are the reconstructed energy and momentum of the B meson candidate 
in the CM frame. We retained events having M bc > 5.2 GeV/c 2 and \ AE\ < 0.2 GeV. The 
signal region was denned as 5.272 GeV/c 2 < M bc < 5.285 GeV/c 2 and \AE\ < 0.025 GeV. 

To suppress the large combinatorial background dominated by the two-jet-like e + e~ — > qq 
continuum process, variables characterizing the event topology were used. We required the 
ratio of the second to zeroth Fox- Wolfram moments [13] R2 < 0.3 and the thrust value of 
the event T < 0.8. Simulation showed that this selection retained more than 95% of BB 
events and rejected about 55% of cc events and 65% of uu, dd and ss events. 



B° ^ D~D+ decays 

The AE and M bc distributions for B° — > D~D + decays obtained after applying all 
selection criteria described above are shown in Fig. 1. 

The beam-constrained mass distribution is well-described using a Gaussian function for 
the signal and an empirically parameterized background shape introduced by the ARGUS 
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Collaboration [14]. The energy difference is described using two Gaussians with the same 
mean for the signal and a linear function for the background. The normalizations, positions 
and widths of the signal are free parameters of the binned likelihood fit. The solid lines in 
Fig. 1 show the result of fits in the interval -0.08 GeV < AE < 0.10 GeV and 5.20 GeV/c 2 < 
M bc < 5.29 GeV/c 2 . The fit to the AE distribution for events in the M bc signal region gives 
a signal yield of iVdata = 1372 ±42 events. The width of the narrower Gaussian that contains 
about 66% of the signal is 8.6 ± 0.9 MeV. From the fit to the M^ c distribution for events in 
the AE signal region we find 1381 ±45 signal events. Due to the simpler parameterization of 
the energy difference distribution and possible peaking background in the M^ c distribution, 
we choose to evaluate all the results from the fits to the AE distribution for events within 
the Mbc signal region. While the level of the combinatorial background is found to be low, 
we used events in the sidebands of the D s and D meson invariant mass distributions to check 
for possible remaining peaking background as described later. 

The efficiency of the selection criteria depends weakly on the reconstructed decay 
channel of the charmed mesons. Efficiencies as obtained from the simulated event sam- 
ple for different combinations of D and D s meson decay channels are used for calcula- 
tion of the overall efficiency The efficiency includes intermediate branching fractions: 
e(D s D) = ^j^eijBiD^BiDj) = (6.31 ±0.88) x 10~ 4 . Here, e itj represents the ef- 
ficiency for reconstructing the event if the D s meson decays through the i-th mode and 
the D meson through the j-th mode. The largest part of the uncertainty (±13.9%) arises 
due to the uncertainties in the intermediate branching fractions [6, 15]. There is also a 
small contribution due to limited MC statistics. Alternatively, one can express the ef- 
ficiency in terms of that for B(D+ — > 0tt + ), relative to which the branching fractions 
for D~ — > K* Q K~ and D~ — > K$K~ are well-measured [6]. In this case we obtain 
e(D a D) = (3.18 ± 0.25) x 10" 2 • B(D+ -> 0vr + )i3(0 -> K + K~). 

The signal peak includes non-resonant Dj — > K + K~n~ decays resulting in the same final 
state. By fitting the AE distributions of a simulated sample we evaluate the fraction of such 
decays with respect to reconstructed B° — > D~D + decays to be r^Kn — (3.46 ± 0.74)%. 
For the final evaluation of the systematic error we include the uncertainty on the branching 
fraction for the non-resonant D~ — > K + K~tt~ decays [6]. This increases the error on tkk-k, 
which becomes (3.46±1.70)%. Due to their small contribution to the signal, the nonresonant 
component is not included in the branching fraction calculation. 

The AE and M^ c distributions obtained using events in the sidebands of both the D 
and D s mesons are in agreement with the observed background under the 5° ^ D~D + 
signal. No peaking is observed for events in the D mass sideband; however, peaking is 
observed in the D s mass sideband as shown in Fig. 2. A peaking component in the Dj — > 
K*°K~ decay mode is due to the three-body decay B° — > D + K*°K~ , previously reported 
by Belle [16]. Since the K*°K~ invariant mass in these decays also populates the region 
under the D~ peak, we must subtract this contribution from the signal. By fitting the AE 
distribution of the D s sideband and normalizing the yield to the area under the D s peak, 
we determine the fraction of such three-body decays within the B° — > D~D + signal peak to 
be T£>k*k — (2.74 ± 0.57)%. Since the method assumes approximately equal contributions 
of B° — > D + K*°K~ decays in the D s sidebands and in the signal region, we assign an 
additional systematic error to tdk*k- From the K*°K~ invariant mass distribution given 
in [16], we estimate a possible relative change in the yield over the relevant M(K*°K~) 
interval (1.77 GeV/c 2 < M(K* K~) < 2.17 GeV/c 2 ) to be ±40%. With the inclusion of 
this additional uncertainty, the fraction r DK * K is (2.74 ± 1.24)%. The signal peak in Fig. 2 
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FIG. 2: AE (left) and Mb c (right) distribution for the D s sideband. The hatched histogram is 
the contribution of events reconstructed in the D~ — ► (f)ir~ and D~ — > K®K~ decay modes where 
no significant peaks are observed. The full line is the result of the fit used to determine the 
B° D+K*°K- contribution. 

also includes a small fraction of D s signal that populates the D s sideband. This contribution 
was evaluated using the simulated sample and subtracted from the fitted number of events 
in the peak. 

Taking into account the contributions from non-resonant D~ — > K + K~ir~ decays (eval- 
uated using MC simulated events) and from B° — > D + K*°K~ decays, we determine the 
number of signal B° — > D~D + in the sample to be [17] N DaD = (1 — r KK7r — r DK * K )N ds ,t& = 
1287 ± 40, where the error is statistical only. 



B° — > D+D~ decays 

The selection optimized for the reconstruction of D~ decays in the B° — > D~D + events 
was also used to search for B° — > DfD~ decays. In this case two D s meson candidates were 
searched for in each event using the decay modes listed above. The M bc and AE variables 
are then calculated. 

An additional source of background in this decay channel is a cross-feed from B° — > 
D~D + decays, where the D + decays into a K~it + it + or K°ir + final state, and one of the 
pions is misidentified as a kaon. Although these events are shifted from zero in the AE 
distribution, they make observation of a small possible signal difficult. In order to reduce this 
contribution, we calculated for each candidate Df — > K + K~n ± or K { $K ± decay the invariant 
mass of the final states particles with the mass substituted for the mass. If the 
invariant mass was within 20 MeV/c 2 of the nominal D ± mass the event was rejected. The 
MC efficiencies for the reconstruction of B° — > DfD~ events for different combinations of 
Df decay modes are determined and used for calculation of the overall efficiency e(D s D s ) = 



7 




AE[GeV] MbJGeV/c 2 ] 

FIG. 3: AE (left) and M bc (right) distributions for the ~B° -> D+D- decay mode. 

J2ij=i,3^i,jB(D S i)B(D S j) = (1.63±0.39) x 10~ 4 . Again, the main uncertainty (±23.1%) arises 
from the uncertainty in the intermediate branching fractions [6, 15]. The same efficiency 
expressed in terms of B{D+ -> 0vr+) reads e(D s D s ) = (0.383±0.040) • {B(D+ -> 0vr + )i3(0 -> 
K+K-)) 2 . 

The distributions of AE and M bc for B° — > DfD~ candidates are shown in Fig. 3. No 
significant signal is observed. The distributions are fitted using a parameterization of the 
DjD + mode. However, the signal AE distribution is described by a single Gaussian and 
the position and the width of the signal are fixed to the values obtained from the MC and 
rescaled by corresponding factors evaluated for the D~D + mode. A binned likelihood fit 
gives Nd s d s = 3.2 ± 2.3 signal events in the AE distribution, where the error reflects the 
statistical uncertainty. The corresponding value for a fit to the M^ c distribution is 4.7 ± 2.6 
events. 

The sidebands were used to check for a possible peaking background. No structure is 
observed in any of the M^ C -AE distributions for events in the D s mass sidebands (see 
Fig. 4). If the D mass region is included in the D s sideband, we observe a clear contribution 
from B° — > D~D + as expected. 

Results 

The number of signal B° — > D~D + events Nc a £) is converted into a branching fraction 
using the efficiency e(D s D) and the number of BB events (assuming equal production of 
B°B° and B + B~ pairs): B(B° -> DjD + ) = (7.42 ± 0.23) x 10" 3 . 

In addition to the above statistical error, we considered several sources of systematic 
uncertainty as listed in Table I. The largest contribution arises due to imprecise knowledge 
of the intermediate branching fractions of D s and D mesons (±13.9%). We assigned a 1% 
relative error for each of the charged tracks used in the reconstructed final states due to 
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FIG. 4: Invariant mass of D^r candidates with hatched sidebands (upper left). AE (upper right) 
and Mb c (lower left) distribution for events with both D s masses in the sideband. AE (lower right) 
for Df signal candidate combined with a Df from a sideband which includes the f T mass region. 



the uncertainty in tracking efficiency. The particle identification efficiency has a relative 
uncertainty of 2.5% per charged kaon and 2.0% per charged pion. An additional 1.4% 
uncertainty is due to the limited statistics of the MC sample used to estimate the efficiency. 
The fraction of true signal events in the sample (1 — r KKn — r DK * K ) is known to an accuracy 
of ±2.2%. Systematics arising from the description of the AE distribution are evaluated by 
comparing the known number of reconstructed B° — > D~D + events in the simulated sample 
with the fitted yield and is found to be 2.2%. Finally, the uncertainty on the number of BB 
events (1.1%) is taken into account. The sum in quadrature of the individual contributions 



9 



TABLE I: Sources of systematic uncertainty in B(B° — > D s D + ) and B(B° — > DfD s ) measure- 
ments. 



B°^D-D+ B°^D+D- 



Source 


Rel. syst. error [%] 


Source 


Rel. syst. error [%] 


L> s and D B's 


13.9 


D s B's 


23.1 


1 racking 


6.0 


Irackmg 


6.0 


Kaon ident. 


7.1 


Kaon ident. 


8.5 


Pion ident. 


6.4 


Pion ident. 


5.2 


MC statistics 


1.4 


MC statistics 


5.3 


rKKn, r DK*K 


2.2 


Fixed parameters 


4.1 


Fitting model 


2.2 


Fitting model 


9.5 


N(BB) 


1.1 


N(BB) 


1.1 


Total 


18.3 


Total 


28.4 



gives a systematic error of 18.3%. The measured branching fraction is thus 

B(B° -> D~D + ) = (7.42 ± 0.23(stat.) ± 1.36(syst.)) x 10~ 3 . (1) 

The same result can be expressed by separating the largest source of systematic un- 
certainty, B(D~ — > 07r~) B((f) — > K + K~), in order to allow the result to be rescaled to 
future improved measurements of this quantity. In this case, the error arising from inter- 
mediate B's reduces to 7.7% [18] and -> D~ D+) ■ B(Dj -> 07r~)i3(0 -> = 
[1.47 ± 0.05(stat.) ± 0.21(syst.)] x 10" 4 . 

We observe no statistically significant signal in the B° — > D+D~ decay mode. From 
the fitted number of events N DsDg we derive — > D+Dj) = (7.1 ± 5.1) x 10" 5 , with 

a statistical error only. Applying the Feldman-Cousins [19] upper limit procedure to this 
value gives B(B° -> D+D~) < 1.6 x 10~ 4 at 90% C.L. 

Evaluation of systematic errors follows closely the treatment of the B° — > D~D + decay 
mode. Individual sources contributing to B(B° — > DfDj) are listed in Table I. The error 
arising from fitting was estimated by comparing the fitted and true number of simulated 
signal events. Additional uncertainty arises due to the use of a fixed width and position 
of the Gaussian function in the fit. The fit was repeated with parameters changed by one 
standard deviation; this resulted in a 4.1% change in B. The overall systematic error is 
found to be 28.4%. We inflate the upper limit by this amount to obtain 

B(B° -> D+D-) < 2.0 x 10~ 4 at 90% C.L. (2) 



Conclusions 

We have measured the branching fraction for B° — > D~D + decays. The measured value 
is B(B° -> D-D+) = [7.42 ± 0.23(stat.) ± 1.36(syst.)] x 10~ 3 , which represents a large 
improvement in accuracy as compared to previous measurements. 

This branching fraction can be used to obtain a value for the CKM matrix element |V^,| 
following [5]. Using B(B° -> Djir + ) from [1], we obtain \V ub /V cb \ = (7.4 ± 1.2) x 10" 2 . The 



10 



error includes experimental uncertainties on the measured branching fractions (the largest 
individual contribution is due to the statistical error on B{B° — > Djn + ) measurement, while 
the uncertainty on the B(D+ — > 07r + ) cancels in the ratio), as well as the uncertainty of 
the theoretical prediction. This result is in agreement with measurements from semileptonic 
decays [6]. 

For B° — > D+D~ decays we find no statistically significant signal. We set an upper limit 
of B(B° -> D+D~) < 2.0 x 10~ 4 at 90% C.L. This result is the first limit on this decay 
mode and is below the prediction of [3] . 
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